GENERAL POSITION OF A PROJECTION AND OF ITS IMAGE 
UNDER A FREE UNITARY BROWNIAN MOTION 
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Abstract. Given an orthogonal projection P and a free unitary Brownian 
motion Y = {Yt)t>o in ^ W*-non commutative probability space such that 
Y and P are *-free in Voiculescu's sense, the main result of this paper states 
that P and YtPY^ are in general position at any time t. To this end, we 
study the dynamics of the unitary operator SYtSY* where S = 2P — 1. More 
precisely, we derive a partial differential equation for the Herglotz transform 
of its spectral distribution, say fit- Then, we provide a flow on the interval 
[—1, 1] in such a way that the Herglotz transform of fit composed with this 
flow is governed by both the Herglotz transforms of the initial (t = 0) and the 
stationary (t = oo) distributions. This fact allows to compute the weight that 
fit assigns to 2 = 1 leading to the main result. As a by-product, the weight 
that the spectral distribution of the free Jacobi process assigns to x = 1 follows 
after a normalization. In the last part of the paper, we use combinatorics 
of non crossing partitions in order to analyze the term corresponding to the 
exponential decay e~"' in the expansion of the n-th moment of SYtSY^. 



1. Reminder 

Let t) be a Vl^*-non commutative probability space with unit 1 and adjoint 
operation *: ^ is a von Neumann algebra endowed with a faithful tracial state r. 
In a recent paper, we studied the free Jacobi process ([9]): this is a family of positive 
operators J = {Jt)t>o valued in the compressed non commutative probability space 



t{P) 

where P G ^ is an orthogonal projection. Actually, the operator Jt is defined by 

Jt = PYtQY^P 

where Y = (lt)t>o € ^ is a free unitary Brownian motion, Q S ^ is another 
orthogonal projection such that {P, Q} and {Y, Y*} are free families in Voiculescu's 
sense (see [T], [2] for details). When P ~ Q and t(P) = 1/2, the main result proved 
in [5] asserts that the spectral distribution of Jt in {Ps:/ P,Tp) fits that of 

Y^t + Y,* + 21 



4 

in (^,t). Two proofs leading to this description were written in [9]. One of them 
is based on the following expansion: let S :~ 2P — 1, t{S) = 0,b = YtSYi* then 
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The description of the spectral distribution of Jt then follows from the fact that in 
this case, SYtSYi and ¥21 share the same spectral distribution ([S], Lemma 1). For 
general ranks t{P) = 6' e [0, 1], we already noticed in [5] that 



(1) r[{PY,PYrr] = ^ + ^ + ^ 5: ( -M{SY,sY:n 



^2" \n - k) ' 



However, we inferred there that a description of the spectral distribution of SYtSY^ 
for arbitrary ranks 9 € (0,1), say /z^, similar to that of /ij^^ seems to be quite 
difficult. In this paper, we shall be mostly interested in the weight that /xf assigns 
to z = 1 and this is two-fold. On the one hand, this weight is exactly the one 
that the spectral distribution of Jt assigns to a; = 1. This follows from the fact 
we prove below that the limit as 71 —> 00 of the RHS of ([T]) depends only on this 
weight, together with Lebesgue convergence Theorem applied to the LHS of ^ 
after normalizing by 1/t(P). On the other hand, a result due to J. von Neumann 
([16|. Theorem 13.7. p. 55) shows that for any orthogonal projections Pi,P2 € 
with subspaces Hi and H2 respectively, 

lim t((PiP2)") =t(Pi AP2), 

n— f 00 

where Pi A P2 denotes the orthogonal projection onto Hi O H2- Accordingly, the 
weight of /i( at z = 1 allows to determine whether or not the projections P and 
YfPY* are in general position at any time t > 0, that isQ: 

(2) r(P A (YtPY;)) ^ max(r(P) + T{YtPYt*) -1,0) = max(2r(P) - 1,0). 

In the case of two projections P,Q with equal rank r(P) = t{Q) = 1/2, it was 
recently proved in [6] that P and YtQY^* are in general position at any time t: 

t{P A (YtQY,*)) = 0. 

This result was then used in order to get an insight into the so-called unification 
conjecture for the projections P, Q. In the same spirit, we shall determine the 
weight of /i( assigned to z = 1 and prove that ([2]) holds for any € (0, 1] and 
any time t > 0. To this end, we use free stochastic calculus to derive a recursive 
time-dependent relation for the moments 

r«(i) := r[{SYtSY;r] = f z-^il{dz), 

T being the unit circle. The obtained relation is then transformed into a partial 
differential equation (hereafter pde) for the Herglotz transform of Ht : 



i/''(t,z):=l + 2^r^(t)z", \z\ < 1. 



Compared to iJ^/^, the pde satisfied by H^ for general 9 conies with the perturba- 
tion 



^If two orthogonal projections are free then they satisfy Q ( 1151 . Lemma 2.1) and are in 
general position according to [B]. 
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while keeping the same initial data H'^{0,-) = i/^/^(0, •)■ Using the method of 
characteristics, we shall find a fiow {t, z) i-> ip^{t, z) on IR+ x [—1, 1] such that 

(3) - [HLi^)r = [H<^{t,^'{t,z))r - [H<^io,z)r. 

Here is the Herglotz transform of the spectral distribution fi^ of SUSU* where 
U & £/ is Haar distributed ([IS])- Equivalently, is the weak limit as i — > od of 

1 /2 

/x^ and is a deformation of the Haar distribution fioo on T (since SYtSY^ and Y2t 
are equally distributed when 6 = 1/2). Up to our best knowledge, no description 
of /Lt^ have showed up yet in literature and it is sufficient for our purposes to focus 
only on its weight fj,f^ at z = 1 (which is the residue of at z = 1). Nonetheless, 
we shall supply here a full description of /x^ relying on an explicit expression of . 
However rather then using the analytic machinery for the multiplicative convolution 
of probability distributions on T f|10]). we found more convenient to us to write 
down by taking the limit in ([1]) as t — )- oo and by using the moment generating 
function of PUPU*P in {Ps/P,Tp) ([8]). Standard analysis arguments show that 
/i^ admits an absolutely continuous part whose support consists of two symmetric 
(with respect to the real axis) arcs that join each other at z = ±1 if and only if 
t{S) = 0. As to its discrete part, it consists of the single point z = 1 with weight 
|t(5')| to z = 1. Coming back to the flow, we shall prove that for any fixed t > 
we can find a real number zf S (0, 1) such that ip^{t, zf) = 1. Since the pole of 
H^{t, •) at z = 1 is exactly the weight of at z = 1 and since this weight can be 
recovered using a radial limit along the real interval (0, 1], then it is given again by 
\t{S)\ according to As a result, taking the limit as n cx) in ([1]) yields the 
main result of the paper: 

Theorem 1. For any t > and any rank t{P) = 6 G (0,1), the orthogonal 
projections P and YtPY^* are in general position. 

It is worth noting that this is a kind of inverse problem since we rely on the 
knowledge of iJ^ in order to guess the behavior of H^{t,z). Conversely, our re- 
sult shows that the weight |t(S')| comes originally from the spectral distribution 
of SYtSY^* and remain mysteriously unchanged up to infinity. This propagation 
phenomenon happens for Jt as well and is a by-product of our main result: (see [7], 
p. 130 for a description of its stationary distribution): 

Corollary 1. Both the spectral distributions of PYtPY^P and of PUPU*P in 
(Px/P, {1/9)t) assign the same weight at x = 1: 

imax[(2^^-l),0]. 

The paper is organized as follows. We first supply a full description of jjP^ and 
derive a closed formula for its moments through Jacobi polynomials. The time- 
dependent recursive relation for the moments r^{t),n > 1 of fj.^ comes next and 
is an instance of a general formula derived in [1]. The obtained relation is then 
transformed it into a pde satisfied by whose dynamics are analyzed using the 
method of characteristics, leading to the flow ^p^ . Once we do, we prove that the 
limit as n — > OD of the RHS of ([1]) depends only on the weight oi fif at z ~ 1 and 
the existence of at any time t > 0. Doing so allows to compute this weight, 
proving thereby our main result and its corollary. We close the paper by giving a 
special interest in the term corresponding to the exponential decay e""* in r^(t). 
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When = 1/2, it is given by a Lagucrrc polynomial ([2], [5]). For general 9 ^ (0, 1), 
we shall see that the difference occurs in the value this term takes at t = 0. Using 
combinatorics of non crossing partitions, we show that this value is the n-th even 
moment of the 1/2-fold convolution of the self-adjoint operator ai — 02, where 
ai,a2 € are two free copies of S. Of course, the resulting convolution is not 
necessarily a probability measure for general 9 G (0, 1] while it is obviously so when 
9 = 1/2 (it reduces the spectral distribution of a since ai and —02 are equally 
distributed). Using the i?-transform machinery ([10]), we can see that the Cauchy- 
Stieltjes transform of the 1/2-convolution of ai — a2 is a root of a third degree 
polynomial that one can express using Gauss hypergeometric functions. 
Henceforth, we shall omit the dependence of our notations on 9 for sake of clarity. 



2. The Stationary distribution 

This section is devoted to the Lebesgue decomposition of the spectral distribution 
/loo of SUSU*, where we recall that U G is Haar unitary operator. More 
precisely, we show that fioo splits into an absolutely continuous part and a singular 
discrete one supported in {1} with weight |r(S')|. To proceed, we shall write down 
its Herglotz transform Hoc- 

[ ^^//oo(d^) = l + 2Vr„z" 

where we set 

r„ :=r„(oo) =T((5?75i7T), ^ > 1- 
This may be done using the free multiplicative convolution of the unitary operators 
S and USU* ([13]) whose common spectral distribution is given by 

66i + {l-6)6_^. 

However we found more convenient to us to deduce from (|T]) and from the 
knowledge of the moment generating function of PUPU*P in P^/P ([8]). The 
issue of our computations is 

Lemma 1. Set k = 29 — 1 = t{S), then 



Hooiz) = ./i + 4k2 



(1-Z)2 

in some neighborhood of the origin. The equality extends analytically to the open 
unit disc. 

Proof. Define 

r(P) 

These are the moments of the stationary free Jacobi process associated with the 
parameters X ~ 1,9 G (O7 1] (see [S] for notations, see also [5]). For instance, 
equation (1) p. 108 in [^ leads on the one side to 



n>l ^ 
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On the other side, we get after summing ([T} over n > 1 

1 r 1 (26- l)z' 1 z" 



n>l 



(5) 



261 
1 

20 



1 



1 



1-2 

(26*- l)z 



^/^^ " 1-2 
Using the following identity whose proof is written in [9] : 



e ^ 22' 

ri>l fe=l 
1 Z*^ " 



2?i 

n — k 



rk 



22k 

k>l n>0 



2n + 2k\ z 



22n- 



E 

ri>0 



27i + 2fc\ z 



22k 



2k 



l2n 



{i + Vi~z) , |z|<l, 



and comparing (jj) and we get 



Er / M„ a/1 - 46'(1 - 9)Z , . X ^ 

r„[a(z)]" = ^^ J_ > -1, with a(z) = -—-=-. 

„>i VI - (1 + VI - 

Finally, recall from [5| that a is invertible from the open unit disc onto some 
neighborhood of the origin, where the inverse is given by 

-11 ^ 4z 
a (z) = 



As a result 



Hooiz) = 1 + 2^ r„z" = 



n>l 



[l + zf 
^1- 461(1 -6l)a-i(z) 

/l + z2 + 2z(l- 861(1 -61)) 



1 + 4k2 



(l-z)2 



as desired. But since 



z 1^ 1 + 4k- 



(1-Z)2 

does not take negative values when z belongs to the open unit disc, then the last 
statement of the lemma is clear. □ 

Corollary 2. The Lebesgue decomposition of the spectral measure of SUSU* 
is given by 

I 

Moo = l^l^l + Jl - -r-j— 1{| sinV'|>|K|} 

y sin ip 

Proof. From the previous lemma, iJoo admits a pole at z = 1, therefore fi^o assigns 
a weight at z = 1 given by 

- lim ^/{^~~7)^~f4KFz = \k\. 

As to the remaining parts of fioo, we first discard the values k = 0, 1. Indeed, we 
know that fioo reduces to the Haar distribution on T when k = <^ 6* = 1/2, while 
when K — 1 6 ~ 1 one has 



so that /ioo = <^i- Hence, assume |k| e (0,1). Then /ioo — |k|<5i is absolutely 
continuous with respect to the Haar distribution in T and its density is given by: 



Indeed 



3? 



1 + z 



1 + 4k2 



z e T. 



(1-z)^ 

(1-^2)(1-.) 



1-z ^22 + 2(2k2-1)z + 1 + |k|(1 + 2) 

has a continuous extension on the boundary T, since 

z + 2(2^2 - l)z + 1 

does not take negative values (its root lie on T) and since the denominator does 
not vanish on the closed unit disc. The proposition now follows from the Poisson 
representation of analytic functions in the open unit disc extending continuously to 
T. Finally 



3? 



1 



1 + 



(1 



= Wi 



sin^ 



L{|sinV'|>|K|}- 

□ 



We close this section with the following closed form of the moments r„,n > 1 
showing that these are somehow averages over (0, \ k\) of special polynomials : 

Proposition 1. For any k e (—1, 1] 

Jo 

where P^'^ is the n-th Jacobi polynomial of parameters (1,0) ([U]. p.254)- 
Proof. Using the generalized binomial Theorem (|14|. p. 47), we write 

/, , 4^2^ ^ ^ ^ i-l/2)k \ Ak^z f _ V- (- V2)fc 



(l-z)2 



k>l 



k\ 



[l-zf 



E 

A:>1 



fc! 



(-4.^z)'= E 



(2fc)„ 



n>0 



where for a; € M, (a;)fc = a;(a; + 1) . . . (a; + fc — 1) is the Pochhammer symbol (|14j. 
p. 45). Inverting the order of summation and identifying coefficients of z", it follows 
that 



1 " 



(-l/2)fc (2fc)„_fc, 



2^ fc! (n-fc)! 

k=l ^ ' 



{-Ak^Y, n>l. 



Writing (2fc)„_A 



and since 



r(ri + fc)/r(2fc), fc > 1, using Lcgendrc duplication formula (|11|) 
r(2fc) = 22^-i(fc-l)!(l/2)fc. 



(-1/2). 



2fc- 1 
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(1/2)., 



one gets 



2\k 



(-1/2), in), i^K^) 

(1/2), {n-k)\{k-\)\ k\ 

^ 2fc- 1 (n-A:)!(fc- 1)! kl 

— 1 {n)k+i (-Ac2)'=+i 



-("-1)!E 



2fc + l(n-l-A:)!(fc + l)! kl 



^ (1 - n), (71 + 1), (^2)^+1 

fl \ 

^ 2fc + l (fc + 1)! fc! 

(1 - n), jn + 1), (/.)^fe-t 
"""^^^^^ 2fc + l (2), fc! 



Finally 

(1 - n)k jn + l)fc (A.)^fc+i _ (n + l)fe(^ 

2fc + l (2), fc! (2), fc! 

= ^^^^"i(l-2-^) = -nt"i(l-2.^) 
(2)„-i ?i 

where the second equality follows from [14], p. 255. □ 
Now, we proceed to the study of /ij. 



3. The time-dependent regime 

3.1. Time-dependent recursive relation. This paragraph is devoted to the 
proof via free stochastic calculus of the following result: 

Proposition 2. Let 

s,,{t) := e"*T{{SYtSY;r) = e"V„(0, n > 1, 

then 

n-l 

dtSn{t) = -n^Sj{t)sn~j{t) + K^n^e"*, n>2. 

Proof. : it goes along the same lines of that of Proposition 1 in [9], with minor 
modifications due to the cancellations = 1 rather than = P. For the reader's 
convenience, we write the whole proof and recall first Theorem 3.4 in [1]: 

Theorem 2. Let n> 1 and define 

/2„(ai, . . .,a2n,t) := e'"*^T{aiYta2Y^* . . . a2n-iYta2nYt) 
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where {ai, . . . , a2n} G is -k-free with Y . Set /o(^, t) :— t{A) for any A ^ then 
dtf2n{ai, . . . ,a2n,t) = - ^ /2n-(i-fc)(ai, ■■ - jak, a;+i, . . . , a2n,t)fi-k{ak+i, ■ .. ,ai,t)- 



l<k<l<2n 
l-k=a[2] 



/2n-(i-fc)-i(ai7 • • • : flfc-i, afcOi+i, ai+2, ■ • ■ , a2n, 0/i-fe-i("i"fc+i, afc+2, ■ • ■ , ai-i,t) 



l<k<l<2n 
l-k-l=Q[2] 

Now, we specialize Theorem [5] to Ofc = 5* for all 1 < A; < n so that f2n = Sn 
and consider n > 2 (for n = 1, the result is derived for instance from [T], p. 923). 
Since both indices fc, I in the first (respectively second) sum in Theorem [2] have the 
same (respectively different) parity, therefore k and Z + 1 in the second sum have 
the same parity and so do I and fc + 1. Accordingly, the first sum does not contain 
terms fo{-,t) while the second does: they correspond to indices / = 2n, k = 1 and to 
l = k+l,l<k<2n—l. Since S*^ = 1 and since r is a trace, then the contribution 
of indices fc = 1, Z = 2n is 

2 nt 

K e . 

For I = fc + l,f < k < 2n — 1, we distinguish two cases: the contribution of 
1 < < 2?i - 2 is 

(2n-2)K2e"* 

while that of fc = 271 — 1, Z = 2n is K^e'""^'*. Thus, the whole contribution of the 
indices k = 1,1 = 2n and ofl<fc<2ri— l,Z = fc + lis 

(6) 2nK^e"K 

Next we write I = k + 2s + I for integer positive values of s and distinguish n ~ 2 
and 71 > 3. If ?i = 2 then there is no additional term in the second sum, while 
if 77 > 3 we separate k = 1 and 2 < fc < 277 — 3. By the same properties of a,r 
mentioned above, the contribution of indices fc = l,Z = 2s + 2,l<s<77 — 2is 

(7) (77-2)K2e"*. 

For the remaining values of 2 < A: < 277 — 3, we distinguish even and odd ones: the 
contribution of indices k = 2j, 1 < j < n — 2, 1 = 2j + 2s + 1 is 

Ji-2n-i-l , . 

(8) E - .v" ^"'^f 

3 = 1 s=l 

while for = 2j + 1, 1 < j < 77 — 2, / = 2s + 2j + 2 we distinguish l<s<77 — j — 2 
and s = 77 — j — 1 . When l<s<n — j — 2 we get 

n-2n-i-2 , . 

(9) Yl E ^v-* = ^v* ^"~ ^^^"~ ^ 

3 = 1 s=l 

while for s = 77 — j — 1 we get 

n-2 

(10) ^K^e"* = (77-2)K2e"*. 

3 = 1 

Coming to the first sum, its contribution is the same as in [9], Lemma 1: 

(11) -nJ2sn-k{t)sk{t). 

fc=l 



The proposition is proved after summing ([5]), (O, © and pH]) . □ 

3.2. Dynamics of the Herglotz transform. Here, we transform the time-dependent 
recursive relation into a pde governing Herglotz transform z i~>- z) of [it- Recall 
that 

and that the moments 

are the coefficients of the expansion of •) as an analytic function: 

= l + 2^r„(t)z", |z|<l. 

n>l 

Using Proposition [21 we readily get: 

Proposition 3. The Herglotz transform H satisfies the equation 
(12) d,H+'-d,H^=2K'^^^, iJ(0,z) = i±i, |z|<l. 

Proof. Elementary computations show the sequence (r„(i))„>i satisfies 
dtnit) = -r,{t) + K^, 

n-l 

dtrn(t) = -7irn{t) - n'^rj[t)rn-j{t) + K^n^ , n>2. 

Consequently, 

dtH = 2^atr„(i)z" 

n>l 

= 2^2 ^ n^z" - 2 ^ nr„(t)z" - 2 ^ n ^ rj (t)r„„j(t)z" 

n>l n>2 i = l 

i>l n>i+l 



n>l 








(1- 








(1- 




z(l- 




(1- 


Z)3 



2 

J>1 



□ 



Remark. The equation (jl2p is a non homogeneous Burgers equation. It allows to 
retrieve the expression of H 00 already obtained in the previous section. Indeed, any 
stationary solution of (|12p is a solution of dtH = 0, that is H{t, z) = H{z) solves 
the first-order ordinary differential equation 



After integrating and taking into account H{0) ~ 1, we get 

H\z) = ^^^j^^ + i^'(O) = HU^) 
(1 - zf 

9 



4. Resolution of the Burgers equation 

In the sequel, we prove that the dynamics of the Herglotz transform is completely 
determined by the initial condition H{0, •), the long-time behavior and some 
characteristics curves. To make easier the computations, we first use the Mobius 
transform 

1 + 2 

z^y = 

1 — z 

which realizes a one-to-one map between the open unit disc and the right half-plane 
{Jftz > 0}. Indeed, this transform replaces the fraction in the RHS of (|12p by a 
cubic polynomial. To see this, set 

F{t,y)^H{t,z), y=^. 

1 — z 

Then F satisfies the equation 

(13) dtF+^{y'-l)dyF'^'^yiy'^l), F{0,y)^y. 

Observe that the stationary solution Hoc reads after this variable change 



The major step toward the proof of our main result is the following theorem: 
Theorem 3. Let F be a solution of the nonlinear equation (|13p . Then 
(14) F'\t, y)) = K'(f>^{t, y) + {l- K')y\ ^y > 0, 

with 



4>\t,y) = l 



where 



(6 + Aev^*)2 -4K2a 

2 , /I 2\ 2 I, ,2 \ A (1 ~ ) 2/ /-I 2\ 

K + {I - K )y , b = a + K , A = — — — p-^ 1 - y 

(1 + y/ay 



In the z- configuration this reads 



(15) + +(1--^) 

with 

(l){t,y) + l y + 1 

Before going into the details of the proof, some remarks are in order. 
Remarks. 

• For ^y > 0, the trajectory t t-^ (pit, y) may cease to exist at some blow-up 
time, say Ti,{y). This corresponds in the z configuration to the time when 
the curve crosses the unit circle at the singularity z = 1. 

• Since the Herglotz transform H(t,-) is analytic inside the unit circle, the 
function Fit, •) is analytic in > 0} and the identity (|14p is meaningful 
wherever the trajectory 1 1— > (pit, y) does not leave this region. As the square 
root is well-defined there then knowing 0^ is sufficient to know p. We point 
out that in general the curves (f> cross the imaginary axis but there are 
some cases where the curves remain confined therein as for instance when 
ye]0,l]. 
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• Observe that 0(0, y) ~ y so that (jl4p can be rewritten under the form 

Fit, 0(i, y)) = i?(0, y), F{t, y) ^ F{t, y) - K^y\ 

This means that the function F is constant along the trajectories. 

• Keeping in mind the expressions of H^o arid H(0,-), is easily seen to 
be equivalent to ([3]). 

• The square root of a is well defined on {5Ry > 0}. 

Proof. Let (f> be the solution of the ordinary differential equation (hereafter ODE) 

(16) a*0=i(02-l)F(t,0), 0(O,y)=y. 

Then differentiating the function Fi : {t,y) i-> F [t , (pit , y)) with respect to t yields 

dtF^{t,y) = \''cj){t,y) {4>^{t,y)~l). 
Therefore solving the pde (IT5|) reduces to the study of the two coupled ODEs: 

r 9,0-i(02_i)Fi, 

(17) <^ 9tFi = 1^2^ (02-1), 

[ 0(O,y) = y,Fi(O,y) = y. 

It is clear that p7|) entails 

FidtFi - K^(j)dt^ = 0. 
Hence, integrating with respect to t yields ([Ti]) 

F2(t,y)-KV(i, 2/) = y'. 

which in turn leads to 



(18) 



9^0=^(02-1)7(1^^ 
4>{0,y) = y- 



Now, we shall solve ((T8| for fixed y > which is equivalent to z e (—1, 1). First, ob- 
serve that 0(i, y) = ±1 are stationary solutions of (|18p for any n and by uniqueness 
of the solution, we deduce that if < y < 1 then is global in time ancd 

|0(t,y)| < l,VteR+. 

More precisely, assume for instance that there exists T > and < y < 1 such 
that 0(T, y) = 1. Then and i i— > 1 solve the Cauchy problem corresponding to 
the data 0(T, y) = 1. Necessarily, 0=1 which is in contrast with < y < 1. As 
a matter of facy, 1 1-> 0(t, y) is non increasing and crosses the right half-plane with 
limit — 1 as i — >■ (X). Similar arguments show that 

y > l=^0(t,y) > l,Vt e]0,T*[. 

Indeed, the lifespan T* is finite due to the cubic power of the non linearity in ([T8| . 
Next, we need to compute the indefinite integral 

dx 



(1 - a;2)^(l-K2)y2+K22.2 



2lf y = 1 then 0=1. 



11 



for real positive x (which is equivalent to e (—1,1)). First, we perform the 
variable change u = 1 — G (^oo, 1) to transform the integral to 

du 



i\/ a — bu + i 



where we set 

a = (1 - K^)y^ + 

b = {2k^ + [1 ~ K^)y') 



C ~ . 



Note that b^ — Aac = (1 — and that the roots of a — bu + cu'^ = lie in [1, oo]. 

Next, we perform the variable change 



^/a{l — vu) = a — bu + I 

and we easily get 

2av — b 



du ~ —2a—, — 5 TT^dv. 

(av^ — cY 



As a result 



du ^ — f dv 1 , 

2^/a / = — In 



2a — bu — 2Jli\/ a — bu + cu 



iVa -bu + cu^ J 2av -b ^ 

But M < 1 so that 2a — bu > 2a — b — (1 — k^)?;^ > and 

(2a - 6u)^ - 4a(a -bu + cu^) = (6^ - 4ac)u^ > 0. 

Consequently 



i\/ a — bu + cju? \/a 
and if U{t, y) = I - (t>'^{t, y), then 



du 1 2a — bu — 2-v/a-\/ a — bu + cu? 

— ■ In ; : . 



J_ 2a-bU{t,y)-2^^a-bU[t,v) + cU^{t,y) _ 
%A \U{t,y)\ 



for some A ~ A{y, k). Equivalently, 



2a - bU{t, y) - 2^^ a - bU{t, y) + c[/2(t, y) 



\U{t,y)\ 

where A = e^"^. Writing this equality as 

2a -[b + e\e^^]U{t, y) = 2^/^^/a - bU{t, y) + cU'^{t, y) 
and raising it to the square we get 

(19) ((6 + eXe^'f ~ 4ac) U{t, y) = 4a(5 + eAe^*) - Aab ^ 4eaAe^*, 

e S {—1,1} being the sign of U. From the observation made before the sign of U 
does not change through the time. To find the value of A we check the preceding 
equation for t = 



A^ + 2e(fe ^)A + 62- 4ac = 0. 

l-y^ 
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• Case < 1; This corresponds to e = 1 and the above equation becomes 

2a 

The discriminant of this polynomial 

16a 



(20) X'^ + 2{b~- 2)A + 62 -4ac = 0. 



IGay"^ 
16ay^ 



{a~b{l^y^) + cil-y^)^) 
{b-2c+cy^) 



Therefore the only solution of ([20| which is not singular at y = 1 is 

2a 2y2y^ 



A = -6 + 



-b + 2 



I — 1 — y^ 

2(a-l) 2(V^-1) 



y2 

-b+2K^ + 2(l-K^) ^ ^ ^ 
1 + Va 

2 



1 + Va 

(i-/^^)V 



(1 + v/a^^ + (1 - ^2)^2) 



• Case y^ > 1: reproducing the same computation yields to 



(1 + v/«2 + (1 _ ^2)j^2)2 



Hence in both cases we get 



(1+ V^^2 + (i„^2)y2)2V ^ 

Finally, (HH) yields 

U{t,y) = —2 = l-02(t,y) 

(fo + Aev^*) - 4ac 

where we simply write A instead of eA and hope there is no ambiguity. □ 

Now we discuss some properties related to the monotonicity of the flow map. 

Proposition 4. Let \ < yi < y2 and denote by T* the lifespan of the trajectory 
t ^ (j){t, yi). Then T* < T* and 

l< 4>{t,y^) < <t>{t,y2), Vte[0,T2^[. 

Proof. The monotonicity of the flow is a consequence of the comparison principle. 
Indeed, denote by < i-)- the trajectory associated to yi. Since 0i > 1, then 

(21) I 9t</<2 = ^(0^- 1)^/(1 + ^2^2^ 

\ 02(0) > 01 (0). 
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Assume that the curves of 0i and (1)2 can intersect and let T be the first time of 
intersection. This means that 

yt e [o,r[, Mt) < Mt). and mt) = Mt). 

Thus necessary 02 (T) < (f)'^ (T) , but from the differential inequality (PT|) we observe 
that 02 (r) > 01 (T), which is impossible. The inequality between the lifespans 
follows easily from the blow up criterion. □ 



5. Proof of Theorem [T] and Corollary [T] 

After this long wave of computations, we proceed to the proof of theorem [T] and 
of its corollary. It consists of two lemmas: the first one gives the limit of the RHS in 
(HI as n — > oo. The second one shows that for any fixed time t, the flow y i— >■ y) 
blows up at some real yt- Equivalently, there exists a real zt such that ip{t, zt) = 1. 

Lemma 2. The following assertions hold true: 

• Let (j) & [0, 27r], then 

i|:(„':',)(-«+--«)=-'"W2)-2^f:)- 

• Let fi be a probability distribution on the unit circle T, then 

• Recall the spectral distribution fit of the unitary operator SYtSY^ . Then 

hrn T[{PYtPY:r] - i [20 - 1 + 

Proof. 

• Using the fact that 
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• Identifying /i with its image under the map z i— >■ arg(2:) e (— tt, tt) we readily 
get 



1 f')n' 

T ^ 



The result follows from Stirling formula 



1 2n\ ,. 1 „ 
hm — — — I = lim —j= = 0. 



• Due to the trace property of r, the spectral distributions of SYtSY^ and 
of Yi SYtS coincide so that pn is invariant under z i-^ z. Hence 

z''nt{dz)= / {z^ + z'')^it{dz) 
Jt 

and the desired limit follows from ([!]). 

□ 

Recall from the previous section that is global in time if and only if < y < 1. 
However, these values of y correspond to — 1 < z < while we need to reach z = 1 
along the real interval (0, 1]. When y > 1, wc shall prove the following 

Lemma 3. 

(1) For any t > 0, there exists Zf €]0, 1[ such that ip(t, Zt) = I, where ip is the 
flow defined in Theorem^ 

(2) For any t > 0, 

lim (1- z)H(t,z) = 2\k\. 

Z-¥l,Z<l 

Proof. (1) Let t > 0, then ip{t, zt) = 1 is equivalent to 

1 + Zj 

2/t)l = +00, yt = - — -■ 

I- Zt 

With regard to the expression of (f) in Theorem [3l we seek y — yt such that 

(6-Aev^*)^-4ac = 
where we recaU that a = (1 — + k^, 6 = + a, b"^ — 4ac = (1 — and 

Note by passing that when k ~ one has c = 0, 5 = a = i/'^, A = (y — + 1) so 

that 

(5 _ Aev^*)' - 4ac = ^ ^eJ'* = 1 ^ ze*(i+^)/(i-^) - 1. 

The last equality states that the E-transform of l2t attains the value 1 (see [2])- 
For general values of k S (—1, 1), we are led to 
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Equivalcntly 

(1 + ^)S/ + (1 - K'fv^iv^ - Ife^^' - 2(1 + _ 1)(^2 ^ ^y^at ^ q 

But since 

2 a — 



y 



1 - k2 



then 



a - 1 



' e2^*-2(a+K2) (1+^)2 ^J^e^* = 



1~K 



.2 



which simphfies to 

(1 + ^f{a - K^) + (a - K2)(a - l)2e2V^* - 2(a + + ^f{a - l)e^' = 0. 

Writing (a — 1) = (y^ — l)(Va + 1) we simpHfy further to get 

(1 + ^^)2(a - + (a - - ife^^* - 2{a + K^){a - l)e^' = 0. 

Gathering terms proportional to , we are led to 

K2[(l+VH)2 + (7^_l)2e2V^t^2(a-l)e^*] =a[(l+V^)2 + (V^-l)2e2^*-2(a-l)e^*] 
and finally to 



Va + 1 
V^- 1 



Hence either 

(22) 

or 

(23) 



V^- 1 
/^+ 1, 



ev^*[v^-|«;|] 



V^- 1 
aA- 1 



For fixed t > each equation admits a unique solution a > 1 and to select the 
suitable solution we refer to the special case n = \. Naturally, the dependence of 
the solution with respect to the parameter k should be continuous but as we will 
see this property is violated by the second formula at the endpoint n ~ 1. Indeed, 
coming back to (fT8|) we see that the solution is given by 

Consequently, 

U[t,y) = l-cf'{t,y)^ 
which becomes singular when 
(24) 



jy' - l)e* 
(y2 - l)e* - 2/2' 



y 



y>i 



y'-i 

But the left-hand side of converges to zero when k goes to 1, whereas the 
right-hand side diverges at this point and thus this formula is not valid for this 
limiting case. Coming to the formula ((22)) . I'Hopital rule leads to 



lim 



\fa — K 



i^iii _ . . _ = lim 



y 
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This gives exactly the equation (|24p and thus we retain 
2) Fix t > 0, there exists yt > I such that 

Km (f>(t, y) = +00. 

y^Vt 

We point out from Proposition |4] that the hfespan of the trajectory starting at 
y G]l,yt[ is larger than the lifespan t of the curve r t-^ (f)(T,yt). Keeping in mind 
F{t, y) = H{t, z) where y = (1 + z)/(l — z), we get 

lim(l-.)i7(M) = 2 lim ^ = 2 lim 

Z^l- y y^y~ (j){t,y) 

But formula entails 



lim 



F^t,cj){t,y)) 



2 



y^y- 4>^{t,y) 



K 



whence we deduce 



lim (1 - z)H{t,z) = 2\k\ 



the lemma is proved. □ 

Proof of theorem [7] and corollary [7J It is a general fact that the singular discrete 
part of corresponds exactly to the poles of H{t, •). Moreover, the weight that 
fit assigns to a given pole can be recovered using radial limits, in particular, fit 
assigns the weight |r| to z = 1 ([4]). Theorem [T] then follows after taking the 
limit as n — )■ 00 in ([1]) and using the last assertion in the first Lemma of this 
section. Corollary [T] follows readily after normalizing ([T]) by 1/t{P) = 1/0 and 
taking the limit as n — >■ c«: the LHS of ([T|) now tends to the weight that the 
spectral distribution of Jt in the compressed space, say vt, assigns to x = 1 since 



□ 



6. Analysis of the moments 

Set e = = [t{S)]^ and recah from Proposition [2] that s„(t) = e"*r„(i),n > 1 
satisfy 

n-l 

dtSn{t) = -n ^ Sj[t)sr,-j{t) + en^e"* 

si(t) = K^e^ + {I- k"^). 
Recall also from [S] that when e = then 

where L^n'^ is the n-th Laguerre polynomial ([H]). In this section, we perform an 
analysis of the 'leading term' in s„(t) when e S [0, 1), that is the term corresponding 
to the fastest decay e~"* of r„(i). Indeed, it is easy to see by induction that 

n 
k=l 
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where P„(e,t) is a polynomial in the variable t whose degree is n — 1. To proceed, 
we first observe that compared to the equation satisfied by s„(t) when e = 0, the 
deformation comes with the factor e"*. Hence, the polynomials Pn{e,t),7i > 1 still 
satisfy 



n-l 



dtPn{e,t) = -n^F,(e,t)P„_,(e,0, n>2 
fe=i 

Pi(e,t) = (l-e). 

However, one needs to compute P„(e, 0) in order to determine the polynomials 
Pn{e,t),n > 2. We shall see in the sequel that while P„(0,0) = s„(0) = 1, P„(e,0) 
changes dramatically when e € (0,1). To this end, we make use of the formula 
(P]. Theorem 14.4) 

7reNC(2n) 

Here NC(2n) is the lattice of non crossing partitions of size 2n, K{n) £ NG{2n) 
denotes the Kreweras complement of tt, c^r is the free cumulant of the 2n-tuple 
(a, a, ... , a) associated with tt and mx(Tr) is the mixed moment of the 2n-tuple 
{Yt,Yf*, . . . , Yt,Yf*) ([13], Chapter IX). Now, the polynomial P„ comes without any 
exponential factor, thus we only need to focus exactly on partitions tt G NC(2ri) 
such that their Kreweras complements Kijr) are non-parity alternating, that is 
each block of K{Tr) lies either in {1, 3, . . . , 2n — 1} or in {2, 4, . . . , 2n} (we identify 
K{tt) sa {1, 2,3,..., 2n}). Indeed, the fc-th moment of Yf is given by 

e-fc*/2i.^(i) (fci)^ fc>l, 
k 

so that the polynomial P„ corresponds to rriK^-n) for which there is no cancellation 
between Y and Y* . According to [13] Exercise 9.42 p. 153-154, the partition tt runs 
over the set NCE(2ri,) of non crossing even partitions (each block of tt has even 
number of elements). Besides, since the constant term of 

equals 1 for any fc > 1 then we end up with 

P„(e,0)= c.AS.S,...,S). 

7reNCE(2Ti) 

We can write this sum as 

P„(e,0)= ^ ^|^C7r(ai - a2,ai - 02, . . . ,ai - 02) 

ireNC(2Ti) 

where ai,a2 € si/ are two free copies of S and \-k\ is the number of blocks of tt. 
Indeed, by freeness of ai and 02 and multi-linearity of free cumulants, one has 

cy(ai - 02, ai - 02, • . . ,ai - 02) = cy(ai, . . . ,ai) -f cy(-a2, . . . , -02) 

for any block V G tt, whence the equality follows. This new way of expressing 
P„(e, 0) hints to the even moments of the 1/2-fold free convolution of the spectral 

18 



distribution of ai — a2 ([S])- Note that when 6 = 1/2 <^ e = then 01,02,-02 are 
distributed according to the symmetric Bernouhi distribution 

hence the 1/2-fold free convolution of the spectral distribution of oi — 02 is still the 
symmetric Bernoulli distribution. Accordingly, we retrieve P„(0,0): 

P„(0,0) = I x^"^[6, + S^,]{dx) = I. 

However when 6 1/2 e the situation becomes rather cumbersome: the 
spectral distribution of oi is given by 

06i + (1 - 9)5_i 

while that of —02 is given by 

{l-0)9Si +9S^i. 
Equivalently, the i?-transform of ai reads 



while that of —02 reads 



yi + 4y{y 

R-aAv) ^ 

near y = Q. It follows that the i?-transform of the 1/2-fold free convolution of 
01 — 02 is given by 



and that its /-^-transform is given by 



Kiy) := Riy) + = + My + ^) + + My - ^) + 2 ^ 
y 4y 
Inverting K (in composition sense) leads to the third-degree polynomial 
y"" - hi{z)y^ + h2iz)y - h3{z) = 

where 

W ^ - 1 

hi(z) 

h2{z) 
h3{z) 



z{z^ - 1) 

+ e - 1 
4z2(22 _ 1) 

1 

4z(z2 - 1)' 



Performing the variables change 

3y=[3h2{z)~{hi{z)f]^^^u + hi{z) 

we recover the reduced form 

3 2[h,{z)]^-9h,{z)h2{z) + 27 h,{z) 

[3/»2(z)-(/ii(z))2]3/2 
19 



Solutions of this equation may be expressed through Gauss hypergeometric func- 
tions 2F1 as shown in [T^] p. 265-266. 
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